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Abstract 

With the aim of generating new constraints on the OZI suppressed couplings of chiral 
perturbation theory a set of six equations of the Roy and Steiner type for the S- and P- 
waves of the irK scattering amplitudes is derived. The range of validity and the multiplicity 
of the solutions are discussed. Precise numerical solutions are obtained in the range E < 1 
GeV which make use as input, for the first time, of the most accurate experimental data 
available at E > 1 GeV for both irK — > irK and tttt — + KK amplitudes. Our main result is 
the determination of a narrow allowed region for the two S-wave scattering lengths. Present 
experimental data below 1 GeV are found to be in generally poor agreement with our results. 
A set of threshold expansion parameters, as well as sub-threshold parameters are computed. 
For the latter, a matching with the ST/ (3) chiral expansion at NLO is performed. 
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1 Introduction 



Scattering amplitudes of pseudo-Goldstone bosons at low energies probe with a unique sensitiv- 
ity the scalar-source sector of Chiral Perturbation Theory (ChPT) [1, 2]. For instance, recent 
progress in the domain of tttt scattering has provided valuable information on the SU{2) chiral 
limit where the masses of the u, d quarks are set to zero. For this purpose, the tttt Roy equations, 
which have been extensively studied in the past [3, 4, 5], were re-analyzed [6] (in particular, a 
formulation as a boundary value problem was developped) and solved numerically [6, 7] (see 
also [8]). These equations constrain the low-energy 7T7r-scattering amplitude by exploiting si- 
multaneously theoretical requirements and data at higher energies. New data on KI4 decays 
from the E865 experiment [9] could thus be studied with the help of the solutions to the Roy 
equations and a bound on the coupling constant £3 of the SU (2) chiral Lagrangian [10] was 
derived for the first time. Constraints on the SU (2) quark condensate were also obtained along 
similar lines [11]. 

In a parallel way, scattering amplitudes involving both pions and kaons at very low energy 
should allow one to unveil features of the SU(3) chiral vacuum, i.e. that in the limit where 
m u , rrid and m s vanish. The structure of the SU(3) chiral vacuum is worth studying for its 
own sake, since SU(3) ChPT provides relations between many low-energy processes involving 
7r-, K- and 77-mesons. In addition, it is interesting to compare SU(2) and SU(3) chiral limits, 
especially in the scalar sector. A sizable difference between the two limits would indicate that 
sea-quark effects are particularly significant in the case of the strange quark [12, 13]. In previous 
works [14, 15] it was shown that the ratio of the pion's decay constant in SU (2) and SU(3) chiral 
limits could be determined from a sum rule based on ttK scattering amplitudes. The deviation 
of this ratio from 1 would indicate a violation of the large- N c approximation. Let us emphasize 
that the latter is often relied upon to attribute values to some 0(p 4 ) couplings arising in the 
scalar sector of the chiral Lagrangian [2, 16]. Our work is motivated by the desire of determining 
from 7r-fT-scattering experimental data as many chiral couplings as possible (in principle, five 
out of the ten independent 0(p 4 ) couplings of the SU(3) chiral Lagrangian [17, 18]), without 
relying on the large- N c approximation. In this paper, we provide the first step of this analysis, 
by deriving the analogue of Roy equations for the ttK system and solving them numerically. A 
simple matching with the SU (3) expansion is performed while more detailed comparisons with 
SU(2) and SU(3) expansions are left for future work. Further motivation for the study of ttK 
scattering can be found in refs. [19, 20]. 

In the case of tttt scattering, Roy observed [21] that general properties of analyticity, unitarity, 
combined with crossing symmetry, lead to a set of non-linear integral equations that the S- and 
P-partial waves must satisfy. A similar program was carried out by Steiner for ttN scattering [22]. 
Given experimental input at high energies (typically E > 1 GeV), Roy-Steiner (RS) equations 
constrain the low-energy behaviour of partial-wave amplitudes. In the present paper, we derive 
and perform a detailed analysis of a system of RS equations for ttK scattering. In this case, 
s — t crossing relates the ttK — ► ttK and the tttt — > KK amplitudes, leading to six coupled 
equations that involve the four ttK S and P partial- wave amplitudes f^ 2 , fi^-, fo^i fi^ 
and the two tttt — > KK amplitudes , g\. Equations of a similar kind have been considered 
earlier [23, 24, 25]. However, some approximations were invoked in the treatment of these 
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equations and, moreover, no accurate experimental input data were available at that time. Since 
then, high-statistics production experiments have been performed for both irK — > ttK [26, 27] 
and 7T7T — > i^X amplitudes [28, 29]. These experiments provide the necessary input data for the 
RS equations with a level of accuracy comparable to the case of irir scattering. Experimental 
data at lower energies should also be available in the near future: the FOCUS experiment [30] 
has demonstrated the feasibility of determining the irK S-w&ve phase shifts at energies lower 
that 1 GeV from the weak decays of D mesons [31], P-wave phase shifts should be measured 
soon in r decays [32] and finally, direct determinations of combinations of S'-wave scattering 
lengths are expected from planned experiments on kaonic atoms [33]. 

The plan of the paper is as follows. After reviewing the notation, we derive the set of RS 
equations that we intend to solve. The setting is similar to a previous work [25] but we differ in 
the number of subtractions used in the dispersive representations. We aim here at an optimal 
use of the energy region where accurate experimental data are available, while avoiding to rely 
on slowly convergent sum rules. After discussing the domains of validity of such equations, we 
explain our treatment of the available experimental input and of the asymptotic regions. Next, 
we start solving the equations. One first eliminates g[j and g\ and the remaining four equations 
then have a similar structure to the irir Roy equations such that recent results concerning 
the multiplicity of the solutions [34, 35] can be exploited. Finally, we turn to the numerical 
resolution and discuss the resulting constraints on the S-wave scattering lengths. Finally, the 
irK amplitudes near and below threshold are constructed and estimates for the 0(p A ) chiral 
coupling constants obtained from matching with the SU(3) expansion are given. 



2 Derivation of the equations 
2.1 Notation 

Let us recall briefly some standard notation [36]. Firstly, we define from the pion and kaon 
masses 

m± = rriK ± m n , S = m 2 K + m 2 , A = m 2 K — m 2 . (1) 

In this paper, exact isopin symmetry will always be assumed. In the isospin limit, there are 
two independent irK amplitudes F I (s,t), with isospin I = \ and I = |. Making use of s — u 
crossing, the / = \ amplitude can be expressed in terms of the / = | one, 

1 1 3 3 3 

F*(8,t,u) = --Fl(s,t,u) + -Fl(u,t,s) . (2) 

It is convenient to introduce the amplitudes F + and F~ which are, respectively, even and odd 
under s — u crossing. In terms of isospin amplitudes, they are defined as 

F + (s,t,u) = -F*(s,t,u) + |F§(a,t,«) 

F-(s,t,u) = ^F l 2(s,t,u) - ^Fl(s,t,u) . (3) 
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The partial-wave expansion of the irK isospin amplitudes is defined as 

F I (s,t) = 16^(21 + l)P l (z s )f l I (s) • (4) 
l 

where Pi{z) are the standard Legendre polynomials and z s is the cosine of the s-channel scat- 
tering angle 

z s = 1 + — — with A s = (s — m 2 + )(s — m 2 _) . (5) 

In a similar way we can expand F + and F~ , and the corresponding partial- wave projections are 
denoted by and fj^(s). The amplitudes can be projected over the partial waves through 

f '^ = u^rl° dt Pl{Zs)FJ • (6) 

lb7rA s J-\ a /s 

The values of the amplitudes at threshold define the S'-wave scattering lengths, with the following 
conventional normalization 

4 = ^/K) ( 7 ) 

(and similarly for in terms of f^(m\)). 

Under s — t crossing, one generates the 1 = and I = 1 -kit — > KK amplitudes, 

G°(t,s,u) = V6F + (s,t,u) 

G l (t,s,u) = 2F~(s,t,u) . (8) 

The partial-wave expansion of the -kit — > KK amplitudes is conventionally defined as 

G\t, s) = 16nV2 ]T(2Z + l)[q n (t)q K (t)} 1 P^gj (t) , (9) 
I 

where the summation runs over even (odd) values of I for I = (/ = 1) due to Bose symmetry in 
the 7T7T channel. In this expression the momenta q n , qx and the cosine of the t-channel scattering 
angle zt are given by 

^) = \f^ii « = Aq ^ K{t) • do) 

The relations between these partial-wave amplitudes and the S-matrix elements are easily 
worked out 

S((s)] = l + 2i^9(s-m 2 + )f l I (s) 

J ttK — >7TK S 

sl(t)} 4< (^(')^))' +1/ % (t _ 4m |, k . (t) . 
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2.2 Fixed-t based dispersive representation 



To derive RS equations, we assume the validity of the Mandelstam double-spectral representa- 
tion [37] from which one can derive a variety of dispersion relations (DR's) for one variable 4 . 
According to the Froissart bound [40], two subtractions are needed at most for F + and one sub- 
traction for F~ (because s — u can be factored out in the latter case) . More detailed information 
about asymptotic behaviour is provided by Regge phenomenology [41], according to which two 
subtractions are indeed necessary for F + while an unsubtracted DR is expected to converge 
for F~ . However, convergence is rather slow in the latter case since the integrand behaves like 
(s')~ 3 / 2 asymptotically. Therefore, we choose to make use of a once-subtracted DR for F~ in 
order to improve the convergence and reduce the sensitivity to the high-energy domain. 

Fixed-t DR's for F + and F~ , with the number of subtractions as discussed above can be 
written in the following form 



1 f°° 

F +(s t ) = c + (t) + - ds' 
-(*) + -/ ds' 



1 1 
— + 



2s' - 2E - t 



lmF + (s',t) . 



F-(s,t) 



(s'-s)(s'-u) X s 



Im F-(s',t) 



(12) 



These expressions involve two unknown functions of t: c + (t) and c~(t). The basic idea for 
determining these functions is to invoke crossing [21, 22], which can be implemented in various 
ways: for instance, one can use fixed-s or fixed- (s — u) DR's. After trying several possibilities, 
we found that DR's at fixed us provide the largest domain of applicability (these relations, 
sometimes called hyperbolic DR's, were exploited in ref. [25]). We start with a special set of 
hyperbolic DR's (more general hyperbolic DR's will be considered later) in which 



us 



The condition above fixes s and u to be functions of t 



(13) 



s = s A (t) = l^2E-t + y/(t- 4m| )(t - 4mf^) 
u = UA (t) = - (2^ - t - y/(t - 4m£)(t - 4m 2 K 



(14) 



According to Regge theory, the function F + (s A ,t) satisfies a once-subtracted DR which is 
slowly converging. Like in the case of the fixed-t DR for F~ , we choose to improve the conver- 
gence by using a twice-subtracted representation. On the other hand, the function F~(sA,t) is 
expected to satisfy an unsubtracted DR which is well converging. Making use of the fact that 
sa(0) = m+, these DR's can be written in the following way 



F + (s A ,t) = 87nn + a£ + b + t+ - 



ds' 



2s' - 2S + t 2s' - 2S - t 



A s / + s't 



A, 



lmF+(s',t' A ) 



4 For the irK amplitude, the existence of fixed-f DR can be established on more general grounds in a finite 
domain of t [38, 39]. 
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+ 



2 i-c 
57T J4r 



dt' 



F (sA,t) Sirm + a | I 



VEirJiml (t') 2 {t'-t) 

1 1 



S A - U A 



2 2 

mi — ra_ 



1 Z" 00 
- / ds' 

7T J m 2 



t 

+ 7T" 



OO 



L A S / + s't X s 
dt' 



-Im ■ 



ImG°(t',4) 
ImF-( S ',4) 



2vr 7 4m 2 t'(f - t) ^ (t' - Aml){t' - Am 
In these equations, we have used the following notation 



(15) 



A 2 



s' A = s A (t'), 4 = 2S-s'-— , 
together with the relation (s' — s A (t))(s' — u A (t)) = X s / + s't . 



(16) 



These representations involve three subtraction constants: the two scattering lengths a^j~, ag 
and an additional parameter denoted b + . Let us now show that the latter can be computed 
through a rapidly convergent sum rule. We notice first that Oq and b + satisfy slowly convergent 
sum rules, 



8irm + a 

2 T 
m + — m_ 



1 



Im F-(s',t' A ) + 



7T Jvn?, A s 



2vr 



dt' 



Im ■ 



G 1 (t' 



im l f ' yj(f -Aml){t' -4m 2 K ) 
b + = — I pmF+( S ',4) + -L- / -^ImG°(t',4) • 

VT J m % \ s > V67T J4m2 (t'J 2 



(17) 



By combining these two sum rules, we can express the parameter b + as a sum rule which has 
better convergence property: 



^ + 87rm + a 



2 2 

m+ — m_ 



71 Jm 



+ 



1 r°° dt 

7T f 



Im 



G (t', s' A 



2J(t' -4m 2 )(t' -4m 2 K )_ 



(18) 



Why does this sum rule converge more quickly ? In the first integral, the combination F + + F~ 
appears, which is the amplitude for the process tt + K~ — > 7r + X _ . The asymptotic region of 
the integrand corresponds to s — ► oo, u — > 0. The amplitude in this region is controlled by the 
Regge trajectories in the u— channel which is exotic, leading to a fast decrease of the integrand. 
In the second integral, the high-energy tail involves the combination ^G°(t', s') — ^G 1 (t',s') 
for t' — > oo and s' — ► 0. The leading Regge contributions are generated by the if** and if* 
trajectories 



lim Im 

t— »oo, s— >0 



_L G o (t)S) _ l G i (t;S )j = /fr„ W _ /3 x *( s )t^( s ) . (19) 



This difference would vanish if Regge trajectories satisfied exactly the property of exchange de- 
generacy. In nature, this property is not exact but it has long been observed to be approximately 
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fulfilled 5 (see e.g. [41] ), which should lead to a significant suppression of the integrand at high 
energies. Therefore, the two integrals involved in eq. (18) are expected to converge quickly, 
providing a determination of b + with only a mild sensitivity to high energies. 

Combining the two dispersive representations eqs. (12) and (15) for the amplitudes F + and 
F~ , the subtraction functions in eqs. (12) get determined in terms of the two 5"-wave scattering 
lengths and we obtain the following representation for the two amplitudes 



F + (s, t) = %irni + a% + b + t + 



1 



TT 



ds' 



oo 



1 



1 



2s' - 2E + t 



+ 



+ 



-I 

t 2 



ds' 



s' — s s' — u \ s > + s't 
2s' - 2S + t 2s' - 2S - t~ 



Im F + {s',t) 



t z f 
JQtt J 4 



A s / + s't 
dt' 



V6rrJ4ml (t') 2 (t'-t) 



F-(s,t) 



8vrm 



rrij 



Ms 



lmG°(t',s' A ) , 
1 



lmF + (s',t' A ) 



m_ 




lm F~(s',t) 



Im F-(s',t' A ) 



2vr J 4m 2 t'(t' - t) ^ (t' - Aml){t' - Am 2 K ) 



(20) 



where the parameter b + is to be expressed in the terms of the sum rule eq. (18). The domain 
of applicability of this representation is limited by the domain of validity of the fixed— t DR's, 
eq. (12). In sec. 3, we will show that the fixed-t DR's hold for t < Am 2 , which enables us to 
perform the projection of eq. (20) on ttK — ► ttK partial waves. We will also need a representation 
which is valid for t > Am 2 in order to obtain equations for the txtx — > KK partial waves. For 
this purpose, we now consider a family of hyperbolic DR's. 



2.3 Fixed us dispersive representation 

Let us consider a general family of hyperbolic DR's for which 

us = b (21) 

is fixed, b is a parameter with (a priori) arbitrary values and should not be confused with the 
subtraction constant b + introduced in the previous section. We write down a twice-subtracted 
representation for F + {sf ) ,t) and a once-subtracted one for F~(sj,,t), 

1 f°° 

F + (s b ,t) = f+(b)+th + (b) + - ds' 

5 The underlying reason for this property is not understood but could be related to the possibility that the 
large- A^ c limit of QCD is described by a string theory [42, 43]. 



2s' -2T, + t 2s' - 2S - t 



M + s't 



Im F + (s',t' b ) 



6 



+ 



2 rc 
57T JAr 



dt' 



ImG°(t',s' b ) , 



F-(s b ,t) 

Sb ~ Ub 



with the notation 



1 roo 

f-(b) + - 2 ds> 



VEirJiml t' 2 (t'-t) 

1 1 



At + s't At 



Im F-(s',t'' 



t 

+ — 



dt 1 



2vr J Ami - 



Im 



si - 



(22) 



i' b = 2£- S '- J 
At = ( s ') 2 - 2Ss' + 6 . 



46 



(23) 



The representations eqs. (22) are a generalization of the DR's eqs. (15) derived for us = A 2 . They 
involve three unknown functions of b: f + (b), f~(b) and h + (b) (which generalize the subtraction 
constants of eqs. (15) ) The two functions f + (b), f~(b) can be determined by matching eqs. (22) 
with the representations eqs. (20) at the point t = (which lies inside their domain of validity). 
Next, the function h + (b) can be expressed as a rapidly convergent sum rule analogous to eq. (18). 
Putting things together, one finally obtains the following representations involving the two S- 
wave scattering lengths aj, as the only arbitrary constants, 



F + (s b ,t) = 8vrm 

- / ds' 

vr J m i 



at + t- 



'■o 



+ 



m\ — m 
2s' - 2S + t 
At + s't 



ImF + (s',t' b ) - ^^Im[F+( S 'X) - F+(s',0)} 



-—l m [F-(s',t' b )-F-(s',0)] 



2 s' — TP t 

— ImF + (s',0) - —lmF-(s',0) 

A s ' A. s ' 



t r°° dt' 

VT J Ami *' 

F~(s b ,t) 8irm + aQ 



lmG°(t',s' b ) ^ GHt',s' b ) 
V6{t'-t) 
1 



2(4 -<) 



Sb ~ Ub 



+ 



TT 



ds' 



IjmF-tfJj _ J-ImF-( s ', ) 



At + s't 



A, 



- _Im[F-( S ',4)-F-( S ',0)] 



+ 



t r 00 dt 7 

2t7 Ami t'(t' - 







Im 



G\t\s' bi 



(24) 



These representations will allow us to perform projections on the t-channel partial waves for 
t > Ami. 
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2.4 RS equations for //(s) 



RS equations can now be obtained by performing the partial-wave projections of the dispersive 
representations obtained above. Projecting eqs. (20) on the / = 0, 1 ttK — ► irK amplitude we 
get the first four equations, 



Refi(s) = k?(s) 

1 Z" 00 r i I 

+- / ^{^(^tOlm^ ^) + 21f, 1 1 ( S ,Olm(/ 1 1 (0} + 

7T J Ami 

Retf(s) = kf(s) 

{ (*"'(7^)v + ^"'H Im/ ' ?(s ° + l K -M^fhs')} 

+- r dt'{K^s,t')lmg° (t') - KkMlmgKt'ri+dfis) . (25) 

The domain of validity in s of these equations is given by eq. (53) below. In these equations, 
the terms k{(s) contain the contributions associated with the subtraction constants, 

, j, , 1 r \ s f , 7, 87rm + an 3s + to 2 \ 

2 327rs y 2 m+ — ml s - m_ J 

His) = ^{b+ + (-3/ + -) *?P±\ ) . (26) 
1V ; 96vrs y V T ml -m 2 _ J v ; 

The equations involve three kinds of kernels Kfi,(s, s'), K^,(s, t'), and Kf v (s, s') (which appear 
only in the driving terms d{). The kernels Kf v read, for I, I' = 0, 1, 

K (s,s') = - Xs+ ^- s) +L(s,s') 

va , j._ 3(X 8 + 2s( S ' + S )) 3(A s , + 2ss'-2A 2 ) 

bsAoA,' A s 



12s 2 (A s ,+2ss'-2A 2 ) + A 2 

-"■li I s ) s ) — ~ 



2s\ s X s i 

+ 3(A, + 2ss> - 2A 2 )(Ay + 2ss' - 2A 2 ) £( ^ g/) 
A.sA.,' 



with 



L(s, S ') = f 



/ A 2 
log(s' + s - 2E) - log s' - 



s 



(28) 
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Next, the kernels KS,,, K®,, (with I' even) read, 



A s \ _ 
Id*) ~~s¥ 



<>(s,t')= 2 ^( q W K ) l 'y s {^g (l + 
n.(M^^,^{( 1 + ^) 1 o g ( 1 ^)- 2 -i(^) 2 }. (29) 



Finally, the kernels K^,, K\ v (with V odd) read 



+ 



2sf 



s(2s - 2S + t') 

a: 



1 + 



2st' 

a7 



6s# 



(30) 



The analyticity properties of the partial-wave amplitudes ff(s) were established in ref. [44]. 
They can be recovered by considering the various kernels. In particular, the circular cut is 
generated by the kernels K[ l i(s,t l ). 

The terms dj (s) are the so-called driving terms in which the contributions from the partial 
waves with I' > 2 are collected 

d{(s) = -Tds'^l (Kf V {s,s') + 1(1-1)^(8,8')) Im/'OO 

71 Jm + l<>2 I V 6 1 

+ ±(-2I + 5)K%,(s,s')Imf}(s')\ 



~ ^E{^2jKa,Olm^(0 + (-3/ + ^)^Wi(*.Olm»2V+i(0} ■ (31) 



+ 



The kernels i^/(s, s') appear in the driving terms only; the first few which are non- vanishing 
read 



K^(s,s') 
K^(s,s') = 



5A, 



s(X s 



-35(A,) 2 s'(ss' - A 
3s 2 (A s 3 



_ TA^ - A 2 )(( S + ^ + A 2 ) - WS) 
Al3(S ' ' j " 3^? 



(32) 



2.5 RS equations for g®(t), g\(t) 

In order to obtain a closed system of equations we now need two equations yielding the real 
parts of <?o(t) and g\(t) valid for positive values of t. They can be obtained from the family of 
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fixed— us DR's of eqs. (24). Using the relation between the cosine of the t-channel scattering 
angle zt and the parameter 6, 

Zt ~ (2E - i) 2 - 4A2 ' {66) 
the projection is carried out by using 



9o(t) = j^j Q 1 dz t F + (s b ,t) 



! 4^2 f 1 2 F-(a 6 ,t) 
Si W = T7 - / dz t z t . (34) 

167T Jo Sb - Ub 



This yields the following two equations for g|J, g{, 

o _ V3m + f . ta \ f p rff Im</8(f) 

3^6 t r 00 rff T i , 

- — / -^ Im 5i(0 

8 Ww t' yu 7 



1 roc r , 

+ E- / , ds ' G+CtyjImZ+^+tG^^yjIm/jT^') +dg(t) . 
i'=0 ^ ■ /m + 



! 2 x /2m ± Op- t /-°° dt'Im^') 



3(m^-m 2 ) 7T Jiml f (f-t) 
+ - / ds'[G^s')lmfr7{s') + G^{t,s')ImK{s')]+d\{t) . (35) 

The two equations (35) together with the four equations (25) form a complete set of Roy-Steiner 
type equations. The domain of validity of the equations for g$, g\ is given in eq. (54) below. 

The equation for g® involves three kinds of kernels: G^,(t,s'), G^,(t, t'). The kernels 
G^,(t, s') have the following form 

(2s' - 2£ + t) , tB v {t,s') 



G+,(t,s') = Vs(2i' + v{ s ,^ x } t/2 M*s>) - {2s ' ™ + t) Mt,*') 

G 0ll (t,s') = V3(2l' + l)^l 



A., 



(36) 



s' 



where P\> are Legendre polynomials and 



G(x) = arCt ™ Hx) , x= w ^ + t , Q t = (t-4ml)(t-4ml). (37) 

We collect below the expressions for the first few of the terms Ai(t, s'), Bi(t, s') 

A (t,s') = l, B (t,s') = -1 

A 1 (t,s') = l, B 1 (t,s') = l 

A 2 (t,s') = l + ^, B 2{ t, s >) = -( 1+ Q±) (38) 



A7' ^,*;--^-r^7 
a u ,n _ 1 , 10s'* , lOg^lzQt + 6s't) R , _ 1 5Q t , 2Q 2 
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Lastly, we quote a few of the kernels G9, 2l ,(t,t'), Gl 2l , +1 (t,t'), 



G° 04 (t, = ^(t' + t- 4S) [3t'(t' - 4£) - 7t(t - 4S) - 64™^] 



256 



7 ^ <3t'(t' - 4S) - 5t(t - 4S) - 32mlm 2 K ] 



384 



(39) 



In the RS equation for g{, eq. (35), one finds two kinds of kernels G lir (t,s') and G\ v {t,t'). The 
kernels G^ v {t,s') have a structure similar to G^,(t, s') encountered above, 



G- l (t,s') = 4V2(2l' + l){ — 



Fix) 



S + t/2) 2 3A 



l — + C v {t,s')} 



with 



F{x) = -5(G(x) - 1) . 



(40) 



(41) 



G(x) is defined in eq. (37) and both F and G are smooth functions around 0. The pieces Ci>(t, s') 
vanish for V = 0, 1 and, for Z' = 2, 3, read 



2a'i , 2s't [(2s'-2i: + t) 2 -9X sf -Us't] 
G 2 {t,s ) = -^2", 3 (M) = ^3 . 

Finally, we display the first few kernels Gj 2 ^ +1 (t,t') 

Gr3(M') = ^(*' + i-4£) 

G\ 5 (t, t') = ^L( t '+t- 4£)[5t'(i' - 4E) - 9t(t - 4£) - Um%m 2 K ] 

These kernels are seen to be polynomials in t, t' . 

The driving terms, d$(t), d\(t), in eqs. (35) have the following expressions 



(42) 



(43) 



4(t) = £ {- / 2 ds ' [G&fr^hzf+W + tGwMhnfrtf) 

l>=2 77 Jm + 

+ - /" ?fc- 2 (^Olm^-2(0 + G^_ 1 (t,«')Im^,_ 1 (0] } 
dKt) = £ {- / 2 ^'G7,(M')W,7( S ') + - / ^^(MOlm^V-i^')} • (44) 



This completes the derivation of a system of equations of the Roy-Steiner type for irK scattering. 
Let us now discuss the domain of validity of these equations. 
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3 Domains of validity 



It is important to assess precisely the domains of validity of the dispersive representations dis- 
cussed in the preceding section. For this purpose, we will adapt the methods reviewed by 
Holder for the ttN system [45]. The discussion is based on the assumption that the scattering 
amplitudes satisfy the Mandelstam double spectral representation [37], i.e., a spectral represen- 
tation in terms of two variables which involves three spectral functions p st (s',t'), p u t(u',t') and 
p us (u',s'). The boundaries of the support of these spectral functions are shown in fig.l. This 
representation and the expressions for these boundaries are obtained from the consideration of 
box diagrams (see for instance [46]). For the irK amplitude, the st boundary is described by 
the two equations 

(t - Aml)\(s, m 2 K ,Aml) - 16m£(s + 3m\ - 3ml) = 

(t - Wml)X(s, m 2 K ,ml) - Urn^s = (45) 

(the ut boundary is obtained by replacing s by u) while the us boundary is defined by the 
following set of equations 

X(u, m 2 K , 4m^)A(s, m 2 K ,ml) — \§m\m 2 K us + 16m^A 2 (m| : — t) = 

X(u, m|,m^)A(s, m\ , 4m 2 ) — \§m\m\us + lQm 2 K l\ 2 {m 2 K — t) = , (46) 

with 

X(x, y, z) = x 2 + y 2 + z 2 - 2xy - 2xz - 2yz . (47) 

Let us consider first the fixed— t DR's. The spectral functions arising in these DR's must be 
real, which implies that the lines of constant t must not cross the double-spectral boundaries. 
From fig. 1 one sees that this condition confines t in the region, 

-48ml < t < Ami , (48) 

where the lower bound comes from the boundary associated with p us and the upper bound from 
the one associated with p st . The second restriction on the domain of validity arises from the 
fact that the spectral function lmF(s',t) is needed in an unphysical region (except if t = 0) 
and must thus be defined using the partial-wave expansion. The domain of convergence of this 
expansion is the large Lehman ellipse (see for instance [46]). In terms of the cosine of the s- 
channel scattering angle z s >, this ellipse has focal points z s > = ±1 and it is limited by the st 
spectral boundary, 

z s r = 1 + 2s ' Tst ^ . (49) 

The function T s t(s) is obtained by solving eq. (45) which describes the stboundaryiox t as 
a function of s. The point — z™ ax of the ellipse corresponds to another value of t given by 
Tg t (s) = —X s /s — T st (s). For each value of s', the convergence of the partial- wave expansion is 
ensured if — z™ ax < z s i < z" iax , i.e., T' st {s') < t < T st (s'). The us boundary provides another 
similar constraint, but it turns out to be weaker than that obtained from the st boundary. 
The conjunction of the two constraints (reality of the spectral functions and convergence of the 
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partial expansion) leads to the fact that the fixed-i dispersion relation for nK scattering is valid 
in the range 

-23ml < t < 4ml ■ ( 50 ) 

A similar discussion can be carried out for the set of dispersion relations with us fixed, 
us = b. Firstly, the criterion that the hyperbolas us = b do not intersect a spectral function 
boundary yields 

-700ml <b< U20ml (51) 

where the lower bound comes from the st boundary and the upper bound from the us boundary. 
For the hyperbolic DR's, the spectral functions ImF (s',t' b ), ImG 1 (t' , s b ) are also needed in 
unphysical regions (unless b = A 2 ), so that the values of b must be restricted to ensure the con- 
vergence of the partial- wave expansion. Considering the Lehman ellipse related to Im F ± (s', t b ) 
restricts the range to 

-100ml <b< AhOml , (52) 
and no further restriction arises from the Lehman ellipse related to ImG 7 (i', s' b ). 

We can now derive the ranges of validity of the RS equations, which are obtained by pro- 
jecting the DR's over partial waves. Let us start with the fixed-t DR's, the projection over ttK 
partial waves is legitimate provided the range of integration of eq. (6) is included inside the 
range of validity in t of the DR's. One deduces that the RS equations for s-channel partial 
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waves (25) are valid for 

3ml < s < 48"4 • (53) 
In the same way, the projection on tttt — ► KK partial waves is allowed only if the range of 
integration of eq. (34) lies within the range of validity in b of the fixed— its DR's. The last two 
RS equations eq. (35) are thus valid for: 

-15ml < t < 70m l ■ ( 54 ) 
The range of validity in t is significantly larger than that in s. This difference stems from Bose 
symmetry, which applies only to the tttt — > KK channels and implies that only even (odd) partial 
waves appear when the isospin is zero (one). Thus, the t-channel projections can be obtained by 
integrating over the limited range < zt < 1, whereas the projection on s-channel partial waves 
requires integrating over the whole range —1<z s <1. One notes that it is possible to project 
the hyperbolic DR's over s-channel partial waves as well. However, the resulting partial-wave 
equations are valid in the range s < 43m 2 , which is somewhat smaller than the range of validity 
of the partial-wave equations obtained from the fixed— t DR's. 



4 Experimental input 

In the previous sections, we have derived a set of RS equations for the s-channel partial waves 
for I = \i \ and Z = 0,1, and the t-channel partial waves for (1,1) = (0,0) and (1,1), which 
we call "lowest" partial waves from now on. Let us consider these equations in the ranges 
m\ < s < s m and Ami < * < t m . The upper limits of which s m , t m (which will be taken such 
that the equations are valid i.e. s rn < 48ml, t m < 70ml) wm be called matching points. A 
simple examination of the RS equations shows that in order to be able to solve for the lowest 
partial waves below the matching points the following input must be provided: 1) the imaginary 
part of the lowest partial waves for s > s m , t > t m , 2) the imaginary parts of the I > 2 partial 
waves above the thresholds and 3) the phases of g® (t) , g\ (t) in the range Ami <t<tm- We will 
discuss below the experimental status of this input. 

For the s-channel partial waves, we choose the matching point at the border of the range of 
validity: 

s m = 0.935 GeV 2 . (55) 

The reason for this choice is that the experimental data available at present comes from pro- 
duction experiments. One expects the precision to decrease as the energy goes down below 1 
GeV. We will see, for instance, that the / = | S'-wave phase shifts seem rather unreliable below 
1 GeV. In the t— channel the range of validity extends, as we have seen, up to t va \ ~ 1.36 GeV 2 
and one could, in principle, choose the matching point anywhere between the KK threshold 
and t va \. In practice, we choose a value slightly above the KK threshold (see sec. 6.1) 

t m = 1-04 GeV 2 (56) 

For the lowest partial waves above the matching point, and for the higher partial waves, we 
exploit experimental data at intermediate energies 

E < ^ = 2.5 GeV (57) 
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and Regge models for E > yfsz. We aim at determining the lowest partial waves below the 
matching point. For this purpose, an additional information is needed concerning unitarity. We 
will make the usual assumption that elastic unitarity holds exactly below the matching points [6] . 
In other terms, in the ttK channel the possible couplings to ttttK and ttttttK are assumed to be 
negligibly small in the low-energy region. For the 5-wave the validity of elastic unitarity was 
observed experimentally up to the rj'K threshold. In principle, the P-wave can couple to the 
ttttK state but no such coupling has been detected for the K* [53], and potentially important 
two-body channels like K*tt, Kp lie above the matching point. Similarly, in the tttt channel we 
assume that the coupling to 4.7T can be neglected below the KK threshold. 

We discuss now the experimental input used to solve the RS equations, before explaining in 
detail their resolution. 

1=3/2 S-wave 

-8 
-10 
-12 
-14 
-16 

•a? -is 

-20 
-22 
-24 
-26 
-28 

0.6 0.8 1 1.2 1.4 1.6 1.8 2 

E (GeV) 

Figure 2: Experimental data from ref. [27] for the I = I S-wave phase shift 




4.1 -kK — ► tiK data 

Phase shift analyses of the ttK — > irK amplitude have been performed based on high-statistics 
production experiments KN — * KirN by Estabrooks et al. [26] and by Aston et al. [27]. Earlier 
results are much less precise and we will not use them in our analysis. The amplitude n + K + — > 
tt + K + which is purely I = | has been measured by Estabrooks et al. [26]. In practice the I = | 
phase shifts are very small in the range E < 2 GeV except for the S'-wave. This phase shift is 
shown in fig. 2 together with our fit, where a simple parametrization with three parameters is 
used 

ta (^w) = iOT < 58 > 

15 



This parametrization is analogous to the one used in ref [19]. Inelasticity is neglected in this 
channel. 

The amplitude tt + K~ — > ir + K~ which involves the following isospin combination 

F c = pa/2 + l F 3/2 ^ (59 ) 

was measured both in ref. [26] and ref. [27] - the latter experiment has better statistics and 
covers a larger energy range. The amplitude F c can be expanded over partial-waves in the same 
way as eq. (4) and refs. [26, 27] provide the phase $;(s) and the modulus a;(s) of these partial 
waves, 

/f(s) = V27Tla i (s)e i *^ . (60) 




Performing a combined fit of the / = | partial waves [26] and of the parameters ai, &i [26, 27] 
one can separate the two isospin partial waves. The data of Aston et al. for the phases $o and 
$1 and our fits are displayed in figs. 3 and 4 respectively in the range 0.9 < E < 1.5 GeV 
(this energy region plays an important role in our analysis). The fits shown here correspond 
to a parametrization of the partial-wave S-matrices as products of Breit-Wigner S'-matrices, 
allowing for inelasticity in the I = \ amplitude to set in at the r)K threshold. Inelasticity is 
found to remain quite small up to E ~ 1.5 GeV. We also tried different fits based on K-matrix 
parametrizations . 
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Figure 4: Experimental data from ref. [27] for the phase $i and the fit used in the calculations. 



The data of Aston et al. and the fits for both ai and for I = to I = 5 and energy up 
to E = 2.5 GeV are shown in figs. 5, 6 and 7. At energies E > 1.8 GeV, ref. [27] found two 
different solutions A and B for the phase shifts, between which we choose sol A (it was pointed 
out in ref. [19] that solution B violates the unitarity bound). These fits allow us to compute the 
relevant integrals up to E = 2.5 GeV. Above that point, we use a Regge-model parametrization 
discussed in sec. 4.3. 

4.2 nix — > KK input 

For our purposes, a key role is played by the I = and I = 1 tttt — > KK amplitudes, which 
can be determined from ttN — > KKN production experiments in the range t > Am 2 K . We 
will make use of the two high-statistics experiments described in Cohen et al. [28] and Etkin 
et al. [29, 47]. The experiment of Cohen et al. [28] determines the charged amplitude 7r + 7r~ — » 
K + K~, thereby providing results for both (/(j and g\. There are several possible solutions but 
physical requirements select a single one, called solution II b in ref. [28]. Close to the KK 
threshold, the presence of the I = 1 phase allows the authors to accurately determine the I = 
phase. The experiment of Etkin et al. concerns the amplitude 7r + 7r - — > K$K$ which is purely 
1 = 0. Because of the absence of the P-wave in this channel, their determination of the phase of 
g§ close to the threshold (where the D-w&ve phase is very small) is likely to be less reliable than 
that of ref. [28]. Their determination of the magnitude of g[] close to the threshold disagrees 
with that of Cohen et al. and also with earlier experiments [48]. Consequently, we make the 
choice to use the results of Etkin et al. only in the range y/t > 1.2 GeV. 
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Figure 5: Modulus and phase of the S- and P -partial-waves amplitudes from ref. [27] and the 
fits in the region 0.9 < E < 2.5 GeV which are used in the calculations. 



Our input for the phase of 5$ is determined as follows. Below the KK threshold this phase 
is identical to the tttt phase shift because of the elastic unitarity assumption. In the range 
2m 7T < E < 0.8 GeV we use solutions of the tttt Roy equations. Simple parametrizations 
were provided recently in refs. [6, 11]. We use the parametrization of ref. [11] together with 
the scattering lengths corresponding to the "extended" fit, with the central values a[j = 0.228, 
Oq = —0.0382. In the range E > 2mx we perform piecewise-polynomial fits of the data of 
refs. [28, 29] and fixing the threshold value to $[j = 200 ± 15 degrees. This range is an educated 
guess based on considering the data of Cohen et al. as well as tttt data. Finally in the range 
0.8 GeV < E < 2niK we perform a fit to the CERN-Munich data as given by Hyams et al. [49] 
and to the polarized target production data recently analyzed by Kaminski et al. [50]. For the 
modulus of <7q, we have performed piecewise polynomial fits to the data of refs. [28, 29]. The 
data and these fits are shown in fig. 8. 

As far as g\ is concerned, we use the experimental determination of the tttt P-wave phase in 
the range 2m n < E < 2mx obtained from the pion vector form factor measured by CLEO [51]. 
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Figure 6: Same as fig.5 for 
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D- and F '-partial waves. 



This determination is compatible with the results of the analysis of tttt Roy equations and has 
a comparable accuracy. At larger energies, we use the experimental results from Cohen et al. 
for the phase and the magnitude of g\. The whole energy range where data are available can be 
fitted using the following form 



(l + ri$(t))*(l + ri«A:(*))3 
{BW(t,m p ) + (P + p iq 2 K (t))BW(t,m p/ ) + ( 7 + j 1 q 2 K (t))BW{t,m p < / )} (61) 

with 

BW-% my)=m v -t- iVtTy 20 ^ + G «® , G P (t) = Vt . (62) 

2G 7r (mf) V Vt J 

Below the KK threshold, q\(t) vanishes and expression (61) reduces to the Kiihn and Santa- 
maria [52] form used in ref. [51]. We take the values of the parameters j3, 7, m p , m p i, m p n 
determined by CLEO and we fit the parameters C, r%, /3\, 71 to the data above the KK thresh- 
old. The data and the fits are shown in fig. 9. 
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Figure 7: Same as fig. 5 for the G- and H -partial waves. 



The amplitudes with I > 2 play a much less significant role in our analysis and are suppressed 
at low energies. They will be described by simple Breit-Wigner parametrizations associated with 
the resonances /2(1200), 72(1525), P3(1690), /4(2050). Masses and partial decay widths of these 
resonances were taken from the PDG [53]. 



4.3 Asymptotic regions 

As discussed above, we can make use of the partial-wave expansion and experimental data up 
to energies E = y/s~2 = 2.5 GeV for the s- as well as the t-channel. Above that point we use a 
description of the amplitudes based on Regge phenomenology. We will content ourselves with 
very unsophisticated models because this energy region turns out to play a very minor role in 
our analysis. In the regime s' — > oo, t fixed, we use the following expression for the amplitudes 
suggested by dual models a la Veneziano [54, 55, 56, 57] (where exact exchange degeneracy is 
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Figure 8: Inputs for the phase of above the inr threshold and its modulus above the KK 
threshold. The data are from refs. [49, 28, 29}. 
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(64) 



For the parameter A and the Pomeron parameters a, b we adopt values inspired by the discussion 
in App. B.4 of ref. [6] with large errors 

A = 14 ± 5 , cr = 5 ± 2.5 mb , 6 = 8 ± 3 GeV" 2 . (65) 

The intercept and slope parameters of the Regge trajectories are determined from the experi- 
mental spectrum of the p and K* resonances, 



K * = 0.352, a p = 0.475, ai = 0.882 GeV" 1 . 



(66) 



For illustration we compare in fig. 10 the imaginary part of F (s, 0) resulting from our fit to 
the experimental data and the Regge asymptotic form with A = 14 . 

Making use of this, it is easy to evaluate the contributions to the various dispersive integrals 
in the range [s2, do]. In the fixed-t DR's we obtain 



F+(s,t)\ 



»2 



2(A a + st) 

(S2) 2 



as2 bt A 

— — exp— + (2 - a P - ait)T(aP + atf) 
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(ftlS 2 ) 



(67) 
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and 



F-(s,t)\ 
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s 2 (3 — a p — ait)T(a p + ait) 



(68) 



In the same manner we can obtain the asymptotic contributions in the amplitudes described 
through hyperbolic DR's 
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(70) 



To derive these contributions, we have used the following expression for ImG J (t',4) in the 
regime where t 1 — > oo 

q;i6 ( 



4=ImG (t',4)U = JlmGVX)U = fT^F^T* ] 1 + ^(- log(«iO + ^(« X *) i 



vrA 



V6 



(71) 
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Figure 10: Comparison of Im F (s, 0) constructed from experimental data and the Regge asymp- 
totic form eq. (63). 



in which an expansion to first order in the parameter b has been performed. 



5 Initial steps in the resolution 

5.1 Solving for g\ 

We have now all the ingredients to solve the set of RS equations. The first step consists in solving 
eqs. (35) for g\. This problem was discussed a long time ago [58, 59] and we recall the main 
ideas here for completeness. Elastic unitarity implies that the phases $>f of these amplitudes 

gl{t)^e^\g((t)\ (72) 

can be identified with the tttt phase shifts df in the unphysical region t < Am? K according to 
Watson's theorem [60] , and therefore they are known in principle. In the physical region t > Am 2 K 
the phases are determined from experiment as was discussed above. 

On the other hand, the modulus of the i-channel partial waves is not known below the KK 
threshold, and must be determined using the equations (35) satisfied by g$ and g\ which have 
the following simple form 

9 §(*) = Ag«) + -/ im| -^*U, 
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gl(t) = A\(t) + ^ [ 
TT J A 



•°° dt'lmg\{t') 

Ami t' t' — t 



(73) 



In sec. 3, we have shown that these relations can be used up to t ~ 1.4 GeV 2 , which includes the 
whole region inaccessible to experiment where |<7§|, \g\\ are needed. The quantities Af(t) are 
analytic functions with a left-hand cut along the negative t axis and no right-hand cut, as can be 
easily verified using eqs. (35) and the explicit form of the kernels provided in sec. 2. Therefore, 
determining the moduli \gj(t)\ in the range 4m 2 < t < t m from eqs. (73) while the phase is 
known is a standard Muskhelishvili-Omnes problem [61, 62]. The most general solution involves 
arbitrary parameters, the number of which depend on the value of the phase at the matching 
point [61]. We have chosen t m to be slightly larger than 4m 2 K . The I = 1 phase <&\{t m ) is lower 
than tt, which implies that the solution for g\ involves no free parameter. The I = phase, as 
we argued in the previous section, satisfies tt < $g(t m ) < 2tt, such that one free parameter is 
involved in the solution. Let us recall the explicit form of the solutions. One first introduces 
the Omnes function 



exp - 

\ TT 



* <5>{{t')dt r 

4m| t'(t' - t) 



= nj R (t) exp[i<Kf (t)0(t - ±ml)9{t m - t)] (74) 



where Qj R (t) is real. Then, the solutions of eqs. (73) read 
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Notice that the integrands are singular when t' 



n\ R (t')t'(t'-t) 

t m , since the Omnes function behaves as 



(76) 



fif fl (f)~|f-*r. 



4>\ (tm) 



(77) 



but the singularity is integrable. When t t m the integrands diverge but this is compensated 
by the factor of flj (t) multiplying the integrals. It can be shown that the solution satisfies 
automatically the first matching condition (details of the proof are given in App. A) 



lim gf(t m ± e)| so i = g{{t 



m ) | input • 



(78) 



Here, go(t), g\{t) are treated in a somewhat different way from that in ref. [15]. In that work, 
an additional subtraction constant was introduced and the values of the subtraction parameters 
were fixed by imposing that the values of <?i(0), go(0) and ^^(O) be equal to the ChPT pre- 
diction at order p 2 . Now, the behaviour around t = is entirely determined by solving the full 
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set of equations with the appropriate boundary conditions - our constraints are dispersive and 
do not rely on ChPT results. 

At this stage, the formulas (75), (76) for g^t), g\(t) involve three parameters: the two 5- wave 
scattering lengths a^ 2 , a^J 2 that appear in the expressions for Afj(t), A\(t) and an additional 
parameter arj. We will now clarify their role. 

5.2 Matching conditions and uniqueness 

Once <7o(i), g\{t) are expressed according to eqs. (75), (76), the set of four RS equations (25) 
becomes a closed set of equations for the four ttK partial waves ff(s), I = 0, 1, I = |, | . 
The structure of these equations is similar to that of irn Roy equations: the kernels consist of 
a singular Cauchy part and a regular part, and elastic unitarity provides a non-linear relation 
between Reff(s) and Im//(s). The equations must be solved with the boundary condition 
that the solution phase shifts must equate the input phase shifts at the frontier of the region 
of resolution (matching condition). Therefore, we can apply the results derived recently [34, 
35] concerning the number of independent solutions in the vicinity of a given solution. The 
multiplicity index of one solution is determined by the values of the input phase shifts at the 
matching point s = s m (with s m ~ 0.935 GeV 2 ). The experimental phase shifts at s = s m lie in 
the following ranges 

< 5l /2 ( Srn ) < |, | < 5\ /2 {s m ) < 7T, 5 3 /2 (s m ) < 0, <5 1 /2 (s m )<0. (79) 

According to the discussion in ref. [35], the multiplicity index in this situation is m = + 1 — 
1 — 1 = —1, to be compared with m = in the case of irir. This means that our situation 
corresponds to a constrained system: a solution will not exist unless the two S'-wave scattering 
lengths lie on a one dimensional curve. 

In practice, however, the phase shift for the / = |, P-wave is extremely small below 1 
GeV and the experimental input is not precise enough to implement matching conditions in this 
channel in any meaningful way (see fig. 17 below). This leads us to treat the / = | P-wave on the 
same footing as the partial waves with I > 2. For instance, the dispersive representations can be 

1/2 

projected on I = 2 and used to compute the real part of f 2 (s) for s < s m while the contribution 
of Im/2^ 2 (s') for s' < s m in the integrands is negligibly small compared to contributions from 
S- and P- waves; it can be evaluated approximately or even ignored 6 . Dropping one matching 
condition, the effective multiplicity index becomes m = for ttK. The fact that the multiplicity 
index vanishes means that solutions should exist for arbitrary values of the two S- wave scattering 

1/2 3/2 

lengths a , a lying in some two dimensional region, and each solution is unique. 

However, not all solutions are physically acceptable. An acceptable solution must satisfy 
the further requirement that it displays no cusp at the matching point [6] . This condition leads 
to constraints on the subtraction parameters. First, let us consider the t-channel, for which 
we choose the matching point t m to be slightly larger than the KK threshold. As discussed 

6 A second argument to neglect the low-energy contribution of the imaginary part of this partial wave is provided 
by the chiral counting Im/j 2 = 0(Im// >2 ) = 0(p 8 ). 
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in the previous section, the solution for g^t) involves one parameter a - While the equality 
Solsoi = Solinput is automatically guaranteed by eq. (75), the solution I sol exhibits a sharp cusp 
at the matching point in general. Therefore, the no-cusp condition fixes the value of a®. The 
same reasoning can be applied to the ttK partial waves: imposing the no-cusp condition to 
the / = \ S- and P-waves provides two equations which should determine, in principle, the two 

scattering lengths Oq 2 , a^ 2 . In other words, given ideal experimental input data 7 with no errors 
in the ranges s > s m and t > t m , one should be able to fix exactly the two scattering lengths 
by solving the RS equations with the appropriate boundary conditions on the values and the 
derivatives of the phase shifts. Obviously, the actual situation is different from that ideal view: 
the input data are known with errors and only for discrete values of the energy, which introduces 
uncertainties on the boundary conditions and thus on the solutions of the RS equations. This 
point will be addressed in the following section. 



6 Numerical solutions and results 
6.1 Numerical determination of the solutions 

We have described how to solve the RS equations for the irn — > KK partial waves. Assuming 
that the input for s > s m is given as well as the input for I > 2 at all energies, our purpose is to 
determine the three phase shifts 

^ /2 (s), S 1 (s)^S 1 1 /2 (s) S\s) ee 6 3 /2 (s) (80) 

in the range m\ < s < s m , so that the Roy-Steiner equations represented symbolically as 

Re/ a (s) = sin(25 a (s)) = <S> a {5 b ,s} (81) 

are satisfied up to a certain accuracy. We introduce a set of N mesh points m\ < Si < s m (N 
was varied between 16 and 30, the results were very stable) and characterize the accuracy of an 
approximate solution by the measure 

e = max|Rer( Sl )-$ a [^]| . (82) 

i,a 

An exact solution, of course, satisfies e = 0. While it is possible to search directly for minimums 
of e, a more appropriate quantity for minimization algorithms is the chi-square 

N 2 

X 2 = EEl Re / a ^)- $a ^ s «]| 2 ' ( 83 ) 

i=l a =0 

which we have minimized using the MINUIT package [63]. Approximations to the irK phase 
shifts S a (s) are constructed in the form of polynomials or piecewise polynomial parametrizations 

7 The data are assumed to be ideal also in the sense that they ensure the existence of a solution to the 
equations [35]. 
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(we tried several forms) similar to that proposed by Schenk [64]. This is essentially the same 
method as in ref. [6] for the irir Roy equations. The parameters are constrained so that the phase 
shifts are continuous at the matching point and the no-cusp condition applies to 5°(s) and 6 1 (s). 
As discussed in sec. 5.2, these additional conditions fix the values of the two 5-wave scattering 
lengths, which are therefore included as two additional parameters in the minimization of the 
chi-square. 

Let us denote by the number of parameters in the representation of 5 a . Taking n^ ' = 3, 
m 1 ) = 2, = 1 we obtain an approximation to the equation with e ~ 5 • 1CP 3 . Adding one 
more parameter with = 2 makes e go down to e ~ 2- 10~ 3 and with still one more parameter, 
= 3, one obtains e ~ 1 • 10~ 3 . This provides good evidence that the approximations are 
converging to a true solution. Seeking a much higher accuracy would be difficult: all integrals 
must be evaluated with a numerical precision better than e, and the computation of the phase 
shifts involve up to three successive numerical integrations (see eqs. (25), (74), (75), (76)). 




-0.8 1 1 1 1 1 1 1 1 1 

0.4 0.5 0.6 0.7 0.8 0.9 1 1.1 1.2 

s [GeV 2 ] 

Figure 11: Left-hand sides of the RS equations. (25) (lines) compared to the right-hand sides 
(points) after minimization in the range m+ < s < 0.93 GeV 2 . 

The accuracy of the solutions is illustrated in fig. 11. In particular, the figure shows that 
the left- and right-hand sides of the RS equations still agree with a good accuracy well above 
the matching point 8 . This constitutes a consistency condition as discussed in ref. [6]. We have 
checked that its fulfilment is a direct consequence of imposing the no-cusp conditions. At this 
level, there is a notable difference between the irir and the irK RS equations. In the case of nir 
scattering [6], it is found that imposing a single no-cusp condition for the P-wave is sufficient 

8 We are then exceeding the strict domain of applicability of the equations but they are still expected to be 
satisfied approximately. 
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to ensure that the no-cusp condition holds to a good approximation for the S'-waves as well, 
and the consistency conditions are well satisfied. In the irK case, we find that it is necessary to 

1 /2 1 /2 

impose no-cusp conditions for the two phase shifts 5 (s) and 5/ (s). In fact, even after doing 

0/9 

so, we find that a (small) cusp remains for the third phase shift <5 (s). This does not represent 
a serious problem, in practice, because this phase shift is not determined very precisely in the 
vicinity of the matching point. 




0.5 1 1.5 2 0.5 1 1.5 2 

t (GeV 2 ) t (GeV 2 ) 



Figure 12: Comparison of the absolute values of g® and g\ obtained from solving the RS equations 
and the corresponding experimental input from ref. [28]. 



Further consistency conditions ought to be considered in the irir — > KK sector. Here as 
well, one expects that the RS equations should be approximately satisfied above the matching 
point. This point is illustrated in fig. 12 which compares the moduli of <?o(i) and g\{t) computed 
from the RS equations to the experimental input for these quantities. Very good agreement is 
observed for g\{f). In contrast, we find that the agreement for <7o(t) is moderately good. In the 
range t > Am 2 K we have checked that the unitarity bound \S w7T ^ K j^\ < 1 is obeyed. Adopting 
a larger value for the matching point t m improves the input-output agreement for t > t m but 
leads to violation of unitarity for t < t m close to the KK threshold. 

Another consistency check can be performed. In the region where t < 4m^, <7o(i) and g\{t) 
can be obtained not only from eqs. (35) which are based on the fixed— us dispersion relations 
but also from the fixed— t ones which are valid in this domain. Both kinds of DR's agree by 
construction at t = 0, the fact that they should continue to agree for negative values of t is not 
trivial and constitutes a check of consistency of the experimental input and of the RS solutions. 
We show these results for t < in fig. 13. 
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t igure 13: Plot of and g{(t) in the region t < obtained by projecting out two different 
kinds of dispersion relations. The lines are obtained from the fixed— us DR's and the points are 
obtained from the fixed— t DR's. 



6.2 Error evaluations and the S- wave scattering lengths 

Our general procedure for evaluating the errors consists in performing variations of the param- 
eters which enter in the description of the input - making use of the errors on these parameters 
and their covariance matrices as provided by running the MINUIT package [63]. The experi- 
mental errors are assumed to be essentially of statistical origin and the errors at different energy 
points are assumed to be independent. Let us discuss first the case of the I = \ S- and P-waves. 
It is clear that this part of the input plays a crucial role as it controls the boundary conditions 
which determine the two 5"-wave scattering lengths. To begin with, one notes that variations of 
the input in the energy region E > 1.5 GeV has a negligibly small influence, so we will consider 
only the energy region yj s m ^ E ^ 1.5 GeV. We have performed two different kinds of fits 
in order to check the validity of the determination of the phase shifts, their derivatives, and 
the errors obtained from varying the parameters at the matching point E = ^/s m . Firstly, we 
perform "global" fits based on a K-matrix parametrization with six parameters for the S-wave 
and seven parameters for the P-wave. These parameters are determined such as to minimize 
the chi-square in the energy region 0.90 < E < 1.50 GeV. Secondly, we have performed "local" 
fits in which one considers separately a small energy region surrounding the matching point 
0.90 < E < 1.1 GeV and the remaining energy region. In the small region we approximate the 
5-wave phase shift by a quadratic polynomial, 

Sl /2 (s) = a + b(s-s m ) + c(s-s m ) 2 , 0.90 < E < 1.1 GeV (84) 
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while for the P-wave we use a linear approximation after subtracting the tail of the K* resonance 



8 J (s) — arctan — 5 = a + bis — s m ) . (85) 

m K * — s 

The results from these two fits concerning the input at the matching point are shown in table 1. 
One observes that the determinations of the phases at the matching point are in good agreement 
as well as that of the errors. The determinations of the derivative of the P-wave agree while 
those of the derivative of the S-wave are only in marginal agreement. In this case, we consider 
the determination from the global fit to be somewhat more reliable as it has continuity and 
smoothness built in. 





phase 


error 


derivative 


error 


1 = global 


46.5 


0.6 


44.1 


5.8 


I = local 


46.2 


0.6 


56.9 


6.6 


I = 1 global 


155.8 


0.4 


148.0 


2.8 


I = 1 local 


156.2 


0.3 


147.4 


2.9 



Table 1: S- and P-waves inputs at the matching point as determined from two different types 
of fit to the data of Aston et al. [27] (see text). Phases are in degrees and their derivatives in 
degreesxGeV^ 1 . 



We can now derive the constraints on the S-wave scattering lengths which arise upon solving 
the RS equations making use of the available experimental input above the matching point. Let 
us first quote some results concerning the errors. Table 2 shows how the errors affecting the 
various pieces of input propagate to the two S-wave scattering lengths. One can see that the 
two main sources of uncertainty are a) the ttK I = \ S'-wave and b) the tttt — ► KK 1 = 
S'-wave. In contrast, the influence of the partial waves with I > 2 (in which the Regge region 
is also included) is rather modest. Finally, this analysis generates the following results for the 
scattering lengths ag, 

m w a /2 ~ 0.224 ± 0.022, aj 2 ~ (-0.448 ± 0.077) 10" 1 . (86) 

There is a significant correlation between these two quantities, the correlation parameter is 
positive and its value is 

pi a = 0.908. (87) 

2 2 

The one-sigma error ellipse corresponding to the above results for the S-wave scattering lengths 

1/2 3/2 

is represented in fig. 14. Our results are compatible with the band obtained for a , a in 
ref. [25]. We find a much smaller allowed region for the scattering lengths simply because we 
have used considerably better experimental input for the S- and P-waves: in the work of ref. [25] 
no data at all were available for E > 1.1 GeV. Predictions from ChPT at 0(p 4 ) for the S-wave 
scattering lengths were presented in ref. [18]. They are recalled below 

m n a /2 = 0.19 ± 0.02 m^aj 2 = -0.05 ± 0.02 (ref. [18]) (88) 
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,1/2 

Jo 


A 1/2 


,3/2 

Jo 


% 


<?! 


I > 2 


10 2 Aa£ /2 


1.89 


0.28 


0.40 


0.79 


0.05 


0.23 


10 2 Aaf 


0.55 


0.09 


0.39 


0.32 


0.14 


0.11 


lO^af - af ) 


1.35 


0.18 


0.10 


0.55 


0.10 


0.15 



Table 2: Sources of error arising from different parts of the input and the resulting errors in the 
determination of the I = scattering lengths in units of m~ L . 



Within the errors these results appear compatible with those from the RS equations. A 
more refined comparison, however, should take the correlation into account. Computing the 
correlation parameter under the same assumptions as used in ref. [18] for the evaluation of the 
errors one obtains the standard error ellipse shown in fig. 14. One observes that the ChPT 
ellipse is very narrow and does not intersect the corresponding error ellipse resulting from the 
RS equations 9 . If one judges from the size of the 0(p 4 ) corrections as compared to the current 
algebra result, it seems not unreasonable to attribute the remaining discrepancy to 0(p 6 ) effects. 
We quote also our results for the two combinations of scattering lengths proportional to o,q , o,q 

77v («o /2 - % /2 ) - 0.269 ± 0.015, m w {a l J 2 + 2og /2 ) ~ 0.134 ± 0.037 (89) 

which are of interest in connection with the irK atom: the square of the first combination 
is proportional to the inverse lifetime of the atom and the sum of the two combinations is 
proportional to the energy shift of the lowest atomic level [65]. The correlation parameter for 
ag , a,Q is also positive and its value is 

p_+ = 0.925 . (90) 

For comparison, let us mention the results for the combinations proportional to , Oq~ in ChPT, 

m w (aj 2 - aj 2 ) ~ 0.238 ± 0.002, m n (aj 2 + 2a /2 ) ~ 0.097 ± 0.047 [ChPT 0(p 4 )] . (91) 

The uncertainty affecting is remarkably small. This, however, could be an artifact of the 
0(p 4 ) approximation. It remains to investigate how 0(p 6 ) corrections affect this result. 

As discussed above, the fact that the two S-wave scattering lengths are determined indepen- 
dently (to some extent) comes from imposing the no-cusp matching conditions. The difference 
of the two scattering lengths, Oq 2 — a%/ 2 , can be determined in an alternative way from a sum 
rule [66] (see eq. (17)). Using this sum rule one finds 

m^al 12 - a 3 Q /2 ) = 0.251 ± 0.014 [sum rule] . (92) 

In the evaluation, we use our results for the RS solutions in the integration regions s' < s m , 
t' < t m . The propagation of the experimental errors is studied in the same way as explained 

9 This particular shape reflects two features of the scattering lengths osq and in the chiral expansion at order 
p 4 : a) they are essentially uncorrelated (the correlation parameter is p_ + ~ —0.15), b) the error on is very 
small because it involves a single chiral coupling (L5) which is multiplied by m* while involves seven chiral 
parameters which are multiplied by ra\m 2 K . 
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Figure 14: Standard error ellipse for the S-wave scattering lengths obtained from solving the 
RS equations with boundary conditions. The corresponding ellipse in the ChPT calculation at 
0(p 4 ) and the current-algebra result are also plotted. 



above. A rather small error is found, but one must keep in mind that the dependence on the 
asymptotic region is significant here and it is difficult to evaluate the error from this region in a 
very reliable way. The central value arising from the sum rule is smaller than what is obtained 
from the RS solution, but the two results are compatible within their errors. We also note that 
the output of the sum rule is significantly influenced by the values of the scattering lengths used 
as input in the integrand. For this reason, the result obtained here differs from the one quoted 
in ref. [15]. 

Before we present the results for the amplitudes in the threshold region, a few remarks are 
in order concerning the intermediate energy region, that ranges from the threshold up to the 
matching point. Experimental data from production experiments exist below 1 GeV, but one 
has to keep in mind the possibility that systematic errors may have been underestimated in this 
energy region in such experiments (a discussion of the tttt case can be found in ref. [67]). Fig. 15 
shows the I = \ P-wave phase shift from the RS equations compared to experiment. The central 

curve correspond to solving with cl^ 2 , d^ 2 taken at the center of the ellipse fig. 14 while the 
upper and lower curves are obtained by using the points on the ellipse with the maximal and 

1/2 

the minimal values for a respectively. The experimental results are seen to deviate from the 
solutions as the energy decreases from the matching point. In particular, the mass of the K* 
which is predicted from the RS equations is 

m K * = (905 ±2) MeV (93) 
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Figure 15: The I = | P-wave phase shift obtained from solving the RS equations. The central 

1/2 3/2 

curve corresponds to solving with a , a taken at the center of the ellipse fig. 14. The upper 

1 /2 

(lower) curves are obtained by using the points with the maximal (minimal) values for a on 
this ellipse. 

1 /2 

(where mpc* is defined such that 5^ (mx*) = 7r /2) is nearly 10 MeV larger than the mass quoted 
in ref. [27] (m^-*o = 896 ± 0.7 MeV). This discrepancy may appear worrying at first sight. It is 
caused, in part, by isospin breaking which is not taken into account by the RS equations. This 
could generate an uncertainty of a few MeV as to the value of the K* mass that should come 
out from solving the equations 10 . Besides, it cannot be excluded that the mass of the K* may 
not be as accurately known as one might believe. The determinations of the K* + , K*° masses 
used by the PDG are all based on hadronic production experiments. Recently, a measurement 
of the K* + mass based on the r decay mode r — > Ks^v T indicated of shift by 4 — 5 MeV as 
compared to the PDG value [68]. In principle, this method is more reliable because it is free of 
any final state interaction problem, but better statistics are needed to clarify this issue. 

The two S-wave phase shifts predicted by the RS equations are shown in fig. 16. For the 
isospin I = \ the RS solution does not exhibit any of the oscillations appearing in the data 
of ref. [26]. For the isospin / = | phase shift, the experimental data for E < 0.9 GeV lie 
systematically below the RS curve, by 2-3 standard deviations. The RS equations also predict 
the I = § P -wave phase shift, the result is shown in fig. 17. This phase shift displays the unusual 
feature that it is positive at very low energy and changes sign as the energy increases. In the 
region around 1 GeV the results are in qualitative agreement with the experimental data of 

10 For instance, the result depends on the input values for m w and mx for which we used m n = 0.13957 GeV, 
m K = 0.4957 GeV. 
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Figure 16: Same as fig. 15 for the 1=2 S-wave phase shift (curves in the upper half of the 
hgure) and the I = | S-wave phase shift (curves in lower half). 



6.3 Results for threshold and sub-threshold expansion parameters 

The behaviour of amplitudes at very small energies is conveniently characterized by sets of 
expansion parameters, which are particularly useful for making comparisons with chiral expan- 
sions. We consider first the set obtained by performing an expansion around the ttK threshold. 
These parameters are conventionally defined from the partial-wave amplitudes as follows 

A Re/ / (s ) = q 2i ( a l + b j q 2 + 4^ + ..) (94) 

with 

2 2 2 2 4 

s = m\ + + ^ + . . . (95) 

Once a solution of the RS equations is obtained, all the threshold parameters are predicted. The 
two S-wave scattering lengths are determined from the matching conditions, as explained above. 
The other threshold parameters may be obtained from the dispersive representation eq. (20) in 
the form of sum rules. These are obtained by projecting the DR's over the relevant partial wave 
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and then expanding the variable s in powers of q 2 . Divergences may appear in this process 
because derivatives are discontinuous at threshold and it must be specified that the limit is to 
be taken from above. This problem is easily handled by computing some pieces of the integrals 
analytically as explained in ref . [6] . The sum rules are evaluated by using RS solutions below the 
matching points and the fits to the experimental data above. For 1 = 0, 1 we have computed 
the parameters ai, b\ and q in an alternative manner by using our solution for Re//(s) for 
three values of s and solving a linear system of equations. The two methods were in very good 
agreement and the results for the threshold parameters are summarized in table 3. The values 
of the P- waves scattering lengths in ChPT at NLO was given in ref. [18] 

mla] 12 = 0.016 ± 0.003 mlaf 2 = (1.13 ± 0.57) 10~ 3 (ref. [18]) . (96) 
Within the errors, these values are compatible with our corresponding results displayed in table 3. 



ChPT expansions of the amplitude are expected to have best convergence properties in 
unphysical regions away from any threshold singularity. The dispersive representations derived 
in sec. 2 allow us to evaluate the amplitude in such regions. A first domain considered in the 
literature is the neighbourhood of the point s = u, t = 0. The following set of expansion 
parameters are conventionally introduced 

F-(s,t) = u{J2C- J f^} (97) 
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T 1 


I = - 
1 2 


(t l\ (t 3\ 


m\ a[ 
m% a l 2 

7 / 


(0.19±0.01)10- i 
(0.47 ± 0.03)10 6 
(0.23±0.03)10- 4 


(0.65 ±0.44)10~ 3 
(-0.11 ± 0.27)10 4 
(0.91 ±0.30)10- 5 


(0.18±0.01)10- i 
(0.48 ± 0.01)10 4 
(0.14±0.01)10- 4 


m% b 2 
m\ b[ 
ml b J 2 


(0.85 ±0.04)10- 1 
(0.18 ±0.02)10- 2 
(-0.14±0.03)10- 3 


(-0.37±0.03)10- 1 
(-0.92 ±0.17)10- 3 
(-0.96±0.26)10- 4 


0.12 ±0.01 
(0.27±0.01)10- 2 
(-0.48±0.02)10- 4 


7 / 

m^ c\ 


(-0.45±0.04)10~ i 
(0.71 ±0.11)10" 3 


(0.18 ±0.02)10- 1 
(0.51 ±0.09)10" 3 


(-0.62±0.06)10- 1 
(0.20±0.03)10" 3 



Table 3: Results for the threshold expansion parameters (see eq. (94) ) computed from the 
dispersive representations and the RS equations solutions. The third column shows the difference 
of the I = \ and the I = | parameters 



where 

s — u , . 

v = . 98 

4.m K 

It is customary to quote the values of the dimensionless parameters C^j which are related to C^j 
by 

C* = {m« + r+ 2 > C±, C- 3 = (m, + r^ C~ . (99) 

The results for the subthreshold expansion parameters are collected in table 4. The table also 
shows for comparison results from ref. [69] , which used fits to the experimental data of Estabrooks 
et al. [26] combined with earlier sets of data (taking into account the data points below 1 GeV 
as well as above). The authors of ref. [69] observed that the low-energy part of the data of 
Estabrooks et al leads to inconsistencies with a dispersive representation of F~ . The agreement 
with our results is reasonable for the coefficients C~ A . For the coefficients the results are 
compatible within the errors, except for the coefficient C^q, for which we find a somewhat larger 
value. Another point of interest is the Cheng-Dashen point s = u, t = 2ml ■ The value of the 
amplitude F + at this point can be related to the kaon sigma-term [70] (see [71] for a recent 
review). We obtain for this quantity 

F + (m 2 K ,2ml) = 3.90 ± 1.50 . (100) 



6.4 Some implications for the 0(p 4 ) chiral couplings 

In this section we present some results on the 0(p A ) couplings of the SU(3) chiral expansion, 
which are easily derived from the subthreshold parameters obtained above. More detailed com- 
parisons between chiral expansions and dispersive representations of the ttK scattering amplitude 
should be performed, but this is left for future work. The expression for this amplitude in ChPT 
at order p 4 was presented in ref. [17]. More specifically, we will make use of a reformulation of the 
expression of ref. [17] in which Fq is expressed in terms of F w only (and not Fjc) as in ref. [15] (a 
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Table 4: Results for the dimensionless subthreshold expansion parameters defined in eqs. (97) 
and (99). The last two columns indicate the chiral order of the leading tree-level contribution 
to each parameter in 577(3) and SU (2) ChPT respectively. 



factor Jkk(s) is missing in eq. (41) of that reference). From this, it is straightforward to obtain 
the chiral formulas for the subthreshold expansion parameters. We present these in numerical 
form below, in which we use the following values for the masses, the pion decay constant and 
the renormalization scale (i 

m n = 0.13957, m K = 0.4957, m v = 0.5473, F n = 0.0924, p = 0.77 (all in GeV) . (101) 

The dimensionless subthreshold parameters then have the following numerical expressions 
in ChPT at NLO 

C+ = 0.14985 + " * [4Li + L 3 - (4L 4 + L b ) + 2(L 8 + 2L 6 )] 
Cf = 0.45754 + «+ f )m * [-(4L, + L 3 ) + 2L 4 ] + 2 <L 5 

TT TV 



C£ = 0.02554 + 



4 



2m! 

_p4 



5 

4Li + L 2 + -L 3 



C + = 1.67285 + 8m | 4 m ' [4L 2 + L 3 ] (102) 



while the subthreshold parameters ■ read 



C - = 8.42568 + —^£5 



C w = -0.02533 - ^§^L 3 . (103) 

In order to lighten the notations we have denoted the renormalized couplings L\(y, = 0.77) 
simply by Lj. It is now easy to solve for the Lj's making use of the results from table 4, the 
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results for L l5 ...,L 4 are collected in table 5. The errors are obtained, as before, by varying all 
the parameters of the fits to the input data and taking into account the correlations. These 
errors appear to be rather small but they only reflect the uncertainty coming from the input 
data. The dominant source of uncertainty in the determination of the Lj's comes from the 
unknown higher-order terms in the chiral expansion, this uncertainty is expected to be of the 
order of 30-40% . This can be seen from the table which shows the results of some alternative 
determinations based on the KI4 form factors [72, 73, 16] and on irK sum rules [15] 11 . We also 
quote the results that we get for L5 and for L$ + 2Lq which have rather large errors 

10 3 L 5 = 3.19 ± 2.40 10 3 (L 8 ±2L 6 ) = 3.66 ±1.52 . (104) 

The coupling L5 is determined, in principle, from Cq but its contribution turns out to be 
suppressed, as it appears multiplied by a factor m\. In order to determine Lg + 2Lq we used 
Cq and the value L5 ~ 1.4 ■ 10~ 3 derived from Fk/F^. The large uncertainty for L$ + 2Lq 

reflects that affecting the coefficient Cq or, alternatively, the uncertainty in the combination of 

1/2 3/2 

scattering lengths a + 2a . This could improve considerably once experimental results from 
ttK atoms are available. Our result for L4, though affected by a sizeable error, agrees with the 
evaluations [74, 75] based on a dispersive method for constructing scalar form factors [76] and 
is suggestive of a significant violation of the OZI rule in the scalar sector. 





nK Roy-Steiner 


ttK sum-rules 


Kh, 0(p 4 ) 


Kl 4 , 0(p«) 


10 3 


Li 


1.05 ±0.12 


0.84 ±0.15 


0.46 ± 0.24 


0.53 ±0.25 


10 3 


L 2 


1.32 ± 0.03 


1.36 ±0.13 


1.49 ±0.23 


0.71 ±0.27 


10 3 


L3 


-4.53 ±0.14 


-3.65 ± 0.45 


-3.18 ±0.85 


-2.72 ± 1.12 


10 3 


u 


0.53 ±0.39 


0.22 ±0.30 




-0.2 ±0.9 



Table 5: Chiral couplings Ll(fi), u. = 0.77 GeV obtained by matching the dispersive results for 
the subthreshold expansion parameters (see table 4) with their chiral expansion at order p A . 
Also shown are the results from ref. [15] (col. 3) as well as those from ref. [16] in which fits to 
the KI4 form factors were perfomed using chiral expansions at order p 4 (col. 4) as well as p 6 
(col. 5). 



7 Conclusions 



In this paper, we have set up and then solved a system of equations a la Roy and Steiner for 
the S- and P-partial waves of the irK — > irK and the irir — > KK amplitudes. These equations 
are necessary consequences of analyticity and crossing, together with plausible assumptions 
concerning the range of effective applicability of elastic unitarity. In treating these equations, 
the approach advocated recently in ref. [6] was followed, which consists in choosing a matching 
point around 1 GeV and enforcing a set of boundary conditions at this point. As input for this 

11 In that paper, terms of order p 6 were dropped in the dispersive representations and the phase shifts used 
below 1 GeV in the sum rules were not constrained to obey the RS equations. 
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analysis, we have exploited for the first time the high-statistics data which are now available 
from KN — > KttN as well as ttN — > KKN production experiments. 

The main result obtained from solving the RS equations together with the boundary condi- 
tions is the determination of an allowed region for the two S'-wave scattering lengths which is 
shown, as a one-sigma ellipse, in fig. 14. This region is much smaller than the ones resulting from 
older analyses, e.g. ref. [25]; this simply reflects the better accuracy of the experimental input 
data used here. Using this result together with the dispersive representations one can determine 
the ttK scattering amplitude in regions of very low energies inaccessible to experiment. We have 
computed a set of sub-threshold expansion parameters and then matched the result with the 
SU(3) ChPT expansion of the amplitude at NLO [17, 18]. This leads to a determination of the 
Gasser-Leutwyler coupling constants L\, L2, L3, and L4. Comparisons with previous results 
is suggestive of significant 0(p 6 ) effects but certainly not so large as to invalidate the SU(3) 
expansion. The bounds that we have obtained for the S-wave scattering lengths constrain the 
combination 2Lq + L$. 

The value of L4 is of particular interest. Since this low-energy constant violates the Zweig 
rule in the scalar channel, its value is related to the role of sea-quark effects and to the link 
between the SU{2) and SU(S) chiral limits [12, 13]. Moreover, it was recently pointed out that 
the value of L4 can be used to discriminate between different assignments for the scalar-meson 
multiplets [77], such a connection was also illustrated in ref. [74]. The value that we found is 
in agreement with the determination based on the scalar form- factors [75, 74, 78] but disagrees 
with the prediction from the chiral unitarization model [79]. More detailed comparisons with 
ChPT expansions should be performed but this is left for future work. At present, the amplitude 
has been computed at order p in the three-flavour expansion and, more recently, in the two- 
flavour expansion [80] (see also [81]). The latter is expected to have better convergence but it is 
less predictive: let us however remark that the expression of the antisymmetric amplitude F~ 
involves only three SU(2) chiral parameters. 

Another topic of interest in connection with ttK scattering is the problem of localizing 
unambiguously a possible k meson (see ref. [82] for a review of the literature). A naive test 
based on the collision time concept [83] applied to our results for the S-wave phase shift gives 
no indication for a resonance. In principle, our results provide an improved and more complete 
input for an analysis such as performed in ref. [82]. 

Dispersive analyses, of course, cannot replace low energy measurements. Much more strin- 
gent constraints on the S-wave scattering lengths could be derived from the RS equations if 
reliable data were available at low energy. For instance, the analysis could be much improved 
soon, once low-energy data on the P-wave phase shifts are obtained from the r — > Kirv T decay. 
In the long term, S-wave phase shifts could be measured in D — > Kirlvi decays [30]. Finally, 
direct measurements of combinations of S'-wave scattering lengths are planned, based on form- 
ing ttK atoms and measuring their lifetime and the shift of the lowest atomic level [33] (see 
refs. [84, 85] for a discussion of related theoretical issues). 
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A Continuity of g\(t) at t = t m 

In this appendix we prove the validity of eq. (78) 

lim g\ (t m ± e)| sol = g{ (t m )\ inpnt (105) 

for g\{t). We will consider the limit from below, t — ► t~, the other limit can be handled in an 
exactly similar way. Let us start from eq. (76) for g\(t) which expresses the solution in terms of 
the input values for the phase &\(t') and the modulus |g{(i')|. 



< 7l 1 (t)| sol = Ai(t) + / 1 (t) + / 2 (t) (106) 



with 



/ 2 (t) - dt jp^, (107) 

which behaviour when t — > i~ has to be investigated. In a first step, one writes as 

^ AKtQsinc&Kt') | /-^-^^AKOsm^KOl 





\ ftm 




/ ^ 




[" rtm+a 




| / < 



ni^tW - *) i 4m 2 nlflWC? -t)J ' 

where a is a small positive number. When t — > i~ the modulus of fi]; goes like 

. *}(tm) 

~ l*"U^~ (109) 

and thererefore vanishes since &\(t m ) > 0. This implies that the second terms of h,2(t) in 
eqs. (108) also vanish when t — > t~ because the integrals multiplied by 0} remain finite. 



/i(t - O - n ^-J tm _ a dt n\ R (t>mti-t) 



. /'*™ +a , , I oj shi $}(£') 

= ni(t) ;J (_ i4 -^_iLi. (no) 

Assuming that a is small enough we can replace tt\ R (t') by its leading behaviour when t' — ► i m 
(eq. (77) ) and we make the same replacement for f)}(i). Next, we perform the following change 
of variables in the integrals 

t' = (t m -t)v + t m , (111) 
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the limits are then expressed in the following way, 

1 f°° dv 

I 1 (t^t-) = A\(t m )eM^ 1 i(tm))sm^\(t m )- / — 

TT Jq V"il 



(1-V) 

1 dv 

h(t^t-)=gl(t m )sm^(t m )-J o -^—r- (112) 



with 

a = MM . (113) 

TT 

The result on the value of <?j(*)lsol a ^ ^ ne matching point follows from the values of the two 
definite integrals [86] (which are well defined for < a < 1) 

1 f°° dv cos(7ra) 1 f°° dv 1 

-/ — - = ) — J - + i, - — - = — - (114) 

TT Jo v a (l — v) sin(7ra) tt Jo v a (l + v) sin(7ra) 

as this implies 

h(t t m ) = -A\(t m ), I 2 (t t~) = g\{t m ) (115) 

which proves the continuity equation (105) for g\ when t m is approached from below. Similar 
arguments can be used to prove continuity when t m is approached from above. Finally, the proof 
is easily generalized to the case of g[] which involves one more subtraction. 
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